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Resistive behaviors at nonzero temperatures (T > 0) reflecting a quantum vortex-glass (VG) 
transition (the so-called field-tuned superconductor-insulator transition at T — 0) are studied 
based on a quantum Ginzburg-Landau (GL) action for a s-wave pairing case containing mi- 
croscopic details. The ordinary dissipative dynamics of the pair-field is assumed on the basis 
of a consistency between the fluctuation conductance terms excluded from GL approach and 
an observed negative magnetoresistance. It is shown that the VG contribution, G V g, to 2D 
conductance becomes insensitive to T at an apparent VG transition field B* g defined at exper- 
imentally accessible temperatures but depends on the repulsive electron-electron interaction, 
and that, only in the dirty limit with no electron-repulsion, it takes a universal value at low T. 
Available resistivity data near B* g are explained based on our results, and an extension of the 
theory to 3D case is briefly discussed. 
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§1. Introduction 

The issue of quantum superconductor-insulator (FSI) transition induced by an applied field in 
disordered superconductors has become one of long-standing problems in condensed matter physics. 
Most of resistivity dataBBi'Bi' suggestive of an FSI transition have been discussed based on a 
scaling hypothesis postulated by Fisher!^ on resistive behaviors at nonzero temperatures (T > 0) 
accompanying a 2D vortex-glass (VG) transition at T = 0. However, no theoretical calculation 
justifying his scenario on resistive behavior at T > has been reported so far. In fact, recent 
experimental studies!' i have often led even to an argument against the 2D VG ordering, i.e., 
superconducting ordering in B ^ 0, although theoretically the 2D quantum VG phase will exist 
more likely than the 3D VG phased does at T > because the former is formally similar to the 
thermal 3D glass phase due to line defectsSS which certainly exists at T > 0. 

Fisher's scenario is based on the phase-only approximation of a GL action. The phase fluctuation 
in a phase-only model for homogeneously disordered (i.e., nongranular) superconductors has no 
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dissipative dynamics, and hence, a nondissipative dynamics was assumed thereS Further, since such 
a phase-only action is essentially equivalent to the quantum boson action at low energy, possible 
vortex phases were identified with the corresponding ones of the boson system. However, as in the 
description of thermal vortex phase diagram at nonzero temperatures]© a more microscopic GL 
approach in which the pair-field yields a dissipative dynamics should be able to explain physical 
properties systematically even at low temperatures as far nontrivial intermediate phase© 
does not occur as a consequence of a phase-only model valid for nongranular films. For instance, 
the so-called bose insulating stateS'S* will correspond to the quantum vortex liquid regime© with 
insulating GL fluctuation conductivity. 

Actually, there are several reasons why the Fisher's argument on the resistive behavior near a FSI 
transition should be reconsidered theoretically. First, the argument^ on a universal conductance 
value at the VG transition field, B = B vg , and on the scaling behavior at T > of resistance is 
based on consideration and calculations of the fluctuation conductance G S (T = 0, oj —* 0) in the 
disorder-tuned case at@ B = 0, where oj — > means taking the dc limit. However, at a T = 
criticality, this conductance need not be equal to the quantity© G S (T — > 0, u> = 0). It is the 
latter which is measurable in real experiments, while it is the former which is expected to take a 
universal value as a consequence of a quantum continuous transition.© Further, the temperature 
range over which the presence of a T-independent conductance is suggested is broad so 

that the dissipative dynamics, neglected in ref.6, of the pair-field does not seem to be negligible, and 
thus, it is unclear to what extent taking the T — > limit is essential. In addition, the assumption!) 
that the derivation of G S (T = 0, uj — > 0) at B = is applicable to the nonzero field case is not 
justified, because the VG critical properties of response quantities are quite different from those for 
the B = superconducting transition. For example, the 3D VG phase due to point-like quenched 
disorder has no static Meissner response to any disturbance of magnetic field, and hence the VG 
transition is not© accompanied by a divergent diamagnetic susceptibility in contrast to the case 
of the normal-Meissner transition. It means that it is not©© justified at all to apply!!* the scaling 
argument on transport quantities in the normal-Meissner transition to the FSI transition case. 

Previously, one of the present authors has pointed out© the possibility that the intervening 
metallic behavior at low T may imply the vortex flow conductance taking a nearly universal value 
along the 2D melting line B m (T) of the disorder- free vortex solid, if B m (T) is insensitive to T 
in the quantum regime. Such a quantum melting line insensitive to T is satisfied in the dirty 
limit© for a s-wave pairing and neglecting a repulsive interaction between electrons. Further, 
the B m (0) is implicitly assumed in this proposal© to lie above a T = VG transition field B V g. 
As demonstrated elsewhere ,0 however, B V g calculated in a mean field approximation in terms of 
microscopic parameter in the s-wave dirty limit© seems to, contrary to the observed trend J30'© 
increase with increasing disorder measured by (.Eft) -1 due primarily to the r-dependence of the 
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corresponding H c2 (0) = fl^(O) = O.560 o T cO r// 2 oc t , where E-p is the fermi energy, I = vpT the 
mean free path, T c q the mean field transition temperature in zero field and in clean limit, and 4>q the 
flux quantum. This is not physically surprising, because not only the thermal fluctuation but also 
the vortex pinning are expected to be enhanced with increasing the microscopic disorder. Hence, 
an inclusion of the electronic interplay between the repulsive interaction and disorder seems to be 
necessary in order to get both B vg and H C 2(0) decreasing with increasing disorder or decreasing 
the film thickness.© Actually, we have judged© that a situation in which the 2D melting field 
B m (0) lies above B vg will seldom happen and hence that a contribution to the conductivity arising 
from vortex pinning effects may not be negligible in explaining consistently the flat resistance 
curve (the intervening metallic behavior). In ref.ll, the VG fluctuation term G vg (T = 0,u — > 0) 
of conductance was considered as a byproduct of transport properties near the 3D thermal glass 
transition due to line disorder and was argued to be a nonuniversal quantity in general. However, 
there, no microscopic details were taken into account, and the measurable quantity G vg (T —* 0, to = 
0) was not examined .0 

In this paper, we examine the vortex-glass contribution G vg (uj = 0) to the dc conductance for a 
current perpendicular to an applied field at nonzero temperatures based on the microscopic study© 
of quantum GL action for the s-wave pairing case. As demonstrated in ref.19, the coefficient of 
dissipative term has a remarkable T-dependence rather in the vicinity of T = 0, implying that the 
premise in previous works!) that one may start from a bosonic model at T = is not justified 
in the present issue. In §2, G vg (T > 0,lu = 0) in the quantum critical region around and at 
the critical field B vg is examined on a general ground. As in ref.19, the ordinary dirty limit 
neglecting an interplay between an electron-electron interaction and disorder will be called merely 
as the dirty limit, and the dynamics of the pair- field at low frequencies is assumed according 
to ref.19 to remain dissipative even in T — » limit. Then, we find that G vg (B vg ) becomes a 
nonuniversal constant depending on a strength of the electron-repulsion but that, in the dirty 
limit with no electron-repulsion, it unusually becomes a universal constant independent of material 
parameters. In §3, it is pointed out that, in B > B vg , the "bosonic" contribution including G vg , 
arising from the GL action, to the conductance approaches zero in T — > limit and that a negative 
magnetoresistance of a superconducting origini* at lower temperatures is provided by "fermionic" 
fluctuation contributions, such as the Maki-Thompson fluctuation term, excluded from the GL 
description. Based on these results, the resistive behavior near B vg is discussed in §3 and compared 
with existing key resistive data in thin films. In §4 we comment on extensions of the present theory 
to 3D systems in low T limit. 

§2. 2D VG Conductance in the Quantum Critical Regime 

Since an FSI transition usually occurs far from the region near the zero-field superconducting 
transition where the vortex-pair excitations play essential roles, a high B approximation will be 
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invoked in which the pair-field in any static vortex state is described in terms of the lowest Landau 
level (LLL) modes which do not accommodate vortex-pair excitations. If first neglecting the random 
potential terms leading to the vortex pinning effects, the 2D quantum GL action on LLL fluctuations 
\P of the pair-field takes the form© 



S'unp — I d r 



/?E(M0) + 7M)l*«(r)| 2 + ^ f Ai|*(r,«)| 



(2.1) 



where 7, U 4 > 0, (3 = l/k B T, tf(r,u) = £ w * 

tj(r)e a; is a Matsubara frequency for bosons, 
and the fact that the squared gauge-invariant gradient Q 2 = (— iV + 2tt/4>qA) 2 is replaced by the 
factor r^ 2 = 2irB/<po after operating any LLL eigenfunction was used. The mean field H C 2(T)-lme 
is defined as /x(0) = 0. In this paper, we focus on the temperature range defined by T < T™ f , 
where© 

~ 0.15T c0 B/#i(0). (2.2) 

As is explained later, this temperature scale arises from the denominator of diffusion propagators, 
and in T < T™ f , /i(0) and U4 become T-independent on cooling, while the T-dependence of 7 
depends remarkably on the presence of an electron-electron repulsive interaction. 

The random potential terms of GL action were studied in ref.19. At high T (< T c q) and low B 
(< H C 2(0)), they may be represented simply in terms of a single random T c term, i.e., as a local 
potential form, while, in low T and high B case of our interest, they become spatially nonlocal 
reflecting the fact that the only microscopic scale measuring the spatial variations of in high B 
and in 2D is the averaged vortex spacing re. The replicated GL action within LLL arising after 
the random-averaging is of the form© 



^ , (Mo) + 7 M) E lrf*W)l 2 + i^p^ 1 E P {a) (K^) P {a) (-K -u) 



p -"■.«/- k 



E^E/o»(<=V«W"'>(-m> 

a' B " k 



(2.3) 



where a and a' are replica indices, N v is the number of field-induced vortices, and (k, u) is the 
Fourier transform of ^^(r, r)| 2 and expressed by 

p& (k, u) = J2 e^- k2 / 4 ^ * (p - k y /2, ut) ^ (p + k y /2, u 1+ u), (2.4) 

where ^(r) = J2 P ^o(p)^o,p( r ) with LLL eigenfunction uo, p ( r ) in a Landau gauge. Further, the 
length scales and were rescaled, respectively, in the manners, r/rg —* r and [3 1 ^ 2 ^ — > 
The function foo(k 2 ) is positive and a regular function of k 2 (and also of k/|k| when the Fermi 
surface is anisotropic) and, at least in the dirty limit, independent!^ of material parameters such 
as Although its detailed functional form is not known even in the dirty limit, just the property 
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that the wavenumber k in /qo is entirely scaled by rg becomes essential in examining the critical 
conductance. The familiar random T c model corresponds to the specific case in which the /in- 
dependences in foo(k 2 ) are neglected. Although this nonlocality in /oo is safely negligible at high T 
and low B where the field is measured through the ratio £o/ r B with GL coherence length in dirty 



limit £o — V^fV^cO' as already mentioned, it cannot be neglected in high B and low T. As well 
as U4, the bare pinning strength U p can be regarded as being T-independentJll' in T < T™ f . 

Until reaching eq.(2.27), we assume in this section the GL coefficients 7, U4, and U p to be T- 
independent as if the GL action is the expression at T = 0. First, to illustrate properties of VG 
fluctuation, we invoke a systematic loop (or 1/M) expansion and focuS® on its lowest order 
(M = 00) terms by, as a mathematical tool, assuming that the complex scalar pair-field ^ has 
M-flavors. Up to the lowest order in U p , the random-averaged propagator Qq of LLL fluctuation 
in this case satisfies© 

(SoM)- 1 = m(o) +7M + TTflESoM - A ^oM- ( 2 - 5 ) 

( R) 

Here the factor A is a coefficient of a renormalized pinning vertex off-diagonal in the replica 
indices and is given by 

= ^iV- 1 ^e- k2 /2 /oo(k 2 )(1 + (Tvge -kV 2) -2 ) (2 . 6) 

Z7Tr B k 

where cr vg = (V \j2 r Kfir 2 B ) ^o( Ld )- The solution of eq.(2.5) is easily found to have the form 

m 1+<m i M+ , (2 . 7) 

where s(\ui\) = 1 — i V g,o/2 + 7|o»|^o(0)/4, and 

t vg ,o = l- A^(g (0)) 2 (2.8) 

measures a distance from a mean field VG transition point -B V g,o- The quantity <5o(0) is selfconsis- 
tently determined by substituting eq.(2.7) into eq.(2.5) with u = 0. Using it again, i vgj o valid near 
-B vgi o is obtained (see eq.(2.12) below). 

Note that, in eqs.(2.7) and (2.8), the pinning strength U p (or A ) and the Matsubara frequency 
uj appear only as the combinations x p = U p (C?o(0)) 2 /(27rr^) and j^uj = 7 C?o(0)<-<->, respectively, 
and that, far below the quantum-thermal crossover line on the LLL fluctuation, defined byB© 

T * T cr = ir^— 2 (2-9) 

(see sec. 2 in ref.14), U4 appears only as the combination X4 = UaQq^)) / {2 , K 2 r 2 B r ) 2 ) : as seen in a vg . 
As a rule of Feynman diagram analysis in LLL0) this is quite general and independent of the 
approximation used above. In fact, the relation x, 

» « (£o(0)) 2 , which is the same form as that 
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occuring in the thermal 2D LLL case, arises because the pinning strength does not carry nonzero 
frequency, while the relation £4 oc ^o(O) is due to the that the purely dissipative quantum 

fluctuation raises the dimensionality of fluctuation by two and behaves like the thermal 4D LLL 
fluctuation. Based on this general rule, it is easy to generalize the above expressions of the lowest 
order in U p to the case with arbitrary pinning strength. To do this, we first note the definition of 
dynamical VG susce ptibility00> 

Xvg(k;wi jW2 ) = AT- 1 f e ik R < ^ (r ) ^ ( r + R) >< ^ (r + R) 2 (r) > 

= {(3N v )- l e~^' 2 e ^-P')fe< ^fowi)^,^) >< <p { P ' + k y ,u; 2 )<p* ( P + k y ,u 2 ) >, (2.10) 
p,p' 

which appears in the expression of G vg , where the overbar denotes the random average. By this 
definition, the irreducible vertex in a diagrammatic representation of Xvg is found to carry a quantity 



A (R) = C^_ A( ) (2n) 



( R*\ 

in the quantum regime T < T cr and consequently, the quantity i vgi o is given by replacing Aq 
in eq.(2.8) by A^ R \ where A( unknown function of x p and X4. Using, for simplicity, 

eq.(2.5) (with replacement Aq^ — > A^), t vgi o near -B V g,o an d in T — ► limit is independent of the 
details of A^ and expressed in the form 

_ 7r 2 7 r| B - B vgfi 
£v g ,o - —r T 5 (A 14) 

up to 0(B — B vg fi). The i? vgj o-value was examined in details elsewhere© 

In real disordered films, the critical fluctuation accompanying the second order VG transition at 
a critical field in T — > limit is expected to be strong enough to change a temperature dependence 
near -B vg .o- To find how behaviors near i? V g,o an d the -B V g-value itself are affected by the VG 
critical fluctuation, let us consider an effective action for the VG fluctuation field, as in the context 
of spin-glass, expressed as a tensor Q aj0l i (w, a/) with a ^ a'. Here Q a , a ' depends on two frequency 
variables,© reflecting that the pinning vertex carrying A^) does not convey nonzero frequency. 
Although a detailed diagrammatic analysis is needed to construct the action consistently with the 
derivation of eq.(2.12), this procedure can be bypassed here. First, it will be easily seen that, when 
performing a Landau expansion of the effective action with respect to Q a ,a', a Q n -order term is 
accompanied by n frequency-summations. Further, according to the diagrammatic rules mentioned 
above, A^ appears everywhere as the combination x p R ^ = A^ (Qo(0)) 2 , which is a constant at 
B vgi Q and in T — > limit according to eq.(2.8) with Aq R ^ replaced by and a modification to 
be done below. Then, the VG effective action should take the form 



S cS (Q) = J d 2 r 
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+«7(7 (fl) M,7 ( * ) k>l; 4 R) ) )Q*u*^i^2) + E ^(f) m/J IIEE«w + i( w i. w i+i) 

/ m>3 j=l Wj ay 

(2.13) 

where the index j = m + 1 implies jf = 1, c m 's are constants, and the function g satisfies 

( Ft} 

g(0,0;x p ) = 0. First, it will be easily seen by considering this action at T = that any renor- 
malization of the Gaussian mass i vg! o due to the interaction (m > 3) terms is accompanied only by 

(Ft} 

the parameter x p and hence that the resulting renormalized one t vg is given by 

t vg ee t vgfi + c fl = 1 + c fl - f ~ *V B B-B vg ^ (214) 

'- y 4 -°vg,0 

where -B V g — -B vg) o(l ~~ c flKi/( 7r2r B7)) is the fluctuation- corrected VG transition field. As usual, a 
nonzero coefficient Cfi is estimated in terms of a UV-divergent self energy term which, at the one loop 
level in eq.(2.13), occurs from a positive Hartree-like contribution of the m = 4 term. Consequently, 
as expected, c/i becomes positive and thus, B vg decreases with increasing the quantum fluctuation 
strength U&/ (rg7) according to the above-mentioned expression. Next, let us consider the action at 
nonzero T (< T cr ) near B vg where |t vg | -C 1. This situation is formally similar to that in 3D systems 
at nonzero (but low) fields near T c q in which the so-called thermal XY scaling^© T — T c q ~ i? 3 / 4 is 
expected below a crossover field (a line defined from the relation 5t s = Stj in ref.21) corresponding 
to T cr in the present issue. In the case, the field or the vortex density measures the inverse-square 
of the microscopic length re, and hence a scaling relation given above is expected by comparing 
the two lengths, r# and the correlation length of the pair-field ~ (T — T c q)~ Uxy with uxy — 2/3, 
with each other. Similarly, in the present issue near B vg , we have only two variables t vg (B) and 
7^//? oc 7/(/3v / A^)), and hence, the only T-dependent scaling variable is, in the present case, 
{£ vg (B,T 0))V fl) /A or equivalents, X vg = U g {B){P/^) l l zv , where z is the dynamical 
exponent of the quantum 2D VG transition, and £ vg (B,T — > 0) = (t vg (B))~ u with an exponent 
v > is the VG correlation length in T — > limit written in unit of r#. Thus, the scaling relation 

7m T ) • <215) 

measuring the field range of quantum VG critical regime at nonzero T, is expected near B vg , 
and the correlation length £ vg expressed in unit of re should have, at nonzero T, the scaling 
form £ vg (B,T) = (t vg (B))^ u S^(X vg ), where the scaling function S$ has the limiting behaviors; 
Sf(X vg — > 0) — > (X wg ) u , S^(X vg — > +oo) approaches a positive constant. For convenience, we will 
assume S^(X vg ) = (X vg /(c^ + X vg )) u leading to the correlation length 

U(B,T) = L(B) + cJ-^Tj j (2.16) 

( Ft} 

and yielding the correct limiting behaviors, where is a positive constant depending only on x p 
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near B vg . Note that eq.(2.16) implies £, vg {B = B vg ,T) ~ T~ x l z . Further, although a line at finite 
T, on which £ vg apparently diverges, is defined in B < B vg as a consequence of the scaling (2.15), 
the appearance of such a line is not peculiar to the present issue: For instance, although the thermal 
3D XY scaling in clean limit suggests a line at nonzero B obeying this scaling, not a continuous 
transition but rather a first order (vortex lattice melting) transition occurs there. Similarly, the 
presence of such a line in the present 2D dirty case near T = does not contradict the absence!) 
of a true 2D VG transition at T > 0. Rather, eq.(2.16) itself will fail with decreasing B away from 
the quantum VG critical regime, since the absence of the 2D VG transition at T > is usually 
guaranteed by a nonperturbative origin in a vortex solid© 

Now, let us examine the VG fluctuation contribution G vg to the total conductance G. Near or 
below -£f C 2(T)-line at nonzero T, G is generally expressed in the form@@) 

G = G n + G s = G n + 5G S + Gfl + G vg , (2-17) 

where G n and G s are, respectively, the quasiparticle contribution and the superconducting part 
of G, and G s consists of the Aslamasov-Larkin (AL) fluctuation term Gfi + G vg and of other 
fluctuation contributions 5G S excluded from the GL description. The AL contribution is further 
divided into the part G vg due to VG fluctuation and the remaining one Gfi, which in clean limit 
becomes the vortex flow conductivity© deep in the vortex liquid region at nonzero T. However, 
Gfi in B > -B vg must vanishB in T — > limit, although this vanishing may not be detected at 
low enough B in systems with weaker quantum fluctuation. Nonvanishing contribution in T — ► 
limit may arise from G n and© 5G S , although discussing a rigorous T = result is beyond the 
scope of this paper. In this section, we focus primarily on G vg . Other fluctuation term 5G S will be 
commented on in §3. 

To study G vg (and Gfi) under a uniform current, we need the spatially averaged supercurrent 
< jn( r ) >r where is the external frequency. In general, < jn(r) > r will take the form 

2-kI 2 _ — 

< jfi(ri) > P1 = -j—P' 1 Yl < ^CM. I w + n h Qj) (Qi + Q2)^(r2)^+n(ri)| ri =r 2 > ri 

= 2 ^ j nl @ — V C(]u\, \u + 0|) yY<Pi(p 3 uj)ip (p, uj + U)(x + iy) + (p* (p, u)ipi(p, u + tt)(x - iy) ), 

(2.18) 

where Qj = —id/drj + 2irA(rj)/<j)Q, (p\(p,u) a fluctuation field in the next lowest Landau level 
(NLL), and ^ the full pair-field prior to decomposing into the Landau levels. The to and depen- 
dences in C(|a;|,|a;-(-f2|) are negligible when considering G vg (and Gfi). Further, in obtaining the 
second line of eq.(2.18), terms consisting only of higher Landau levels were neglected by assuming 
a high field approximation0'0'0) in which a static vortex state is described within LLL. In partic- 
ular, it is important to note that the NLL modes cannot participate in describing any static vortex 
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solid for a reason of symmetryElP Since this massive NLL modes inevitably appear in considering 
a response to a uniform current, some terms associated with NLL modes have to be added in the 
GL action, which are expressed by 



^; = EE(ftM)"Vl a) (p^)i 



Un 



EE 



Airr%N v , 



^/n(fc 2 )4 a) (k,0) P r*(k,0) 



,09)*/ 



+ (^/ 01 (A: 2 )pS a) (k,0)^)*(k,0) + c.c. 



where k- = k v + ik x , and 



Pi 



(a) 



(k, 0) = £ e^- k2 /VS a) * (P - ^/2, w) # (P + *vA 



(2.19) 



(2.20) 



ui,p 



For convenience of presentation, the NLL propagator £/i(cj) =< \tp^ (p, uj)\ 2 >= (fii(0) + j^\lo\)^ 1 
is assumed here to have been renormalized in a proper way. The bare mass /ii(0) is of order unity 
near i7 c2 (0). 

The terms G/z + G vg of conductance are obtained in terms of Kubo formulaB 

2 



TTTg 
Ml 



Rq{Gji + G vg ) 



<9 \^ iV _^_ 1 y, < y (>,a;)^0 / ,a;)yi0 / ,a; + O)^(p,a; + O) > 



(Si(0)) 2 



(2.21) 



where Rq = Trh/2e 2 is the resistance quantum. In the present high field approximation, Gi(0) 
appears in both G vg and Gfi only as the combination Mi = 2tt1 2 C(0, 0)£7i(0). Below, it will be 
shown that, deep in the vortex liquid regime, this factor takes the universal value 



Mi = vrr j 



(2.22) 



independent of microscopic details as a consequence of gauge-invariance. Deep in the vortex liquid 
regime, the AL conductance Gfi in the pinning- free case can be assumed to be given by the vortex 
flow expression.© Since the vortex flow in a pinning-free system is not affected by the vortex- 
solidification at T m (B), we can focus on the mean field result of the vortex flow conductance, which 
is most easily derived in terms of a harmonic action!!) 

2\ 



_ < |^MF(r)| 2 >, 

O^har - T-2 



EItMIslP + CSiCo))- 1 

n 



2tt 



sl + —(5 Ax z)r% 



(2.23) 



where sl denotes a uniform displacement of vortices, and $mf is the mean field solution of vortex 
solid. The combination r^ 2 (5xsi) + (27r/0o)5A in eq.(2.23) is a consequence of the gauge-invariance 
and ensures the Josephson relation E = —dSA/dt = — x B, where vx = dsi/dt is the vortex 
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velocity. By substituting Bsl = z x 5 A minimizing the above second term into the first term, 
the vortex flow conductivity is obtained as the coefficient of thus obtained first term, which is 
clearly independent of Gi(0). On the other hand, using eq.(2.21) (with no pinning disorder effect) 
in the mean field approximation where (po 's cj-dependence is absent, one obtains a conductance 
expression oc Mj 2 . Under the condition that these two expressions are identical with each other, 
eq.(2.22) follows. In fact, using the general fact,© already used in writing down eq.(2.23), that the 
uniform fluctuation belonging to NLL around the mean field solution ^mf is nothing but a uniform 
displacement solution with amplitude sl )V — isj^, the averaged supercurrent (2.18) is found to be 
proportional to 2irl 2 C(0, 0)(B x + 5 A) up to 0(sl) and 0(5 A). By identifying this with that 
following from eq.(2.23), we again obtain eq.(2.22). Although, due to a pinning disorder effect on 
the static NLL mode, the pinning-free result (2.22) may be subject to a subtle change, we argue in 
terms of eq.(2.8) that, as far as Gi(0) <C Go(0) is safely satisfied, eq.(2.22) is quantitatively valid 
near B vg in the present high field case. 

Before examining G vg further, a form of Xvg needs to be determined. In the Gaussian (or mean 
field) approximation, it is written as 



t , 0<gvg,o) 2 , 2(C vg ,o) 4 7 (fl) M | 2(C vg , ) 4 7 (fi) k 2 | , (2 21) 



Xv g ,o(k;wi,u;2) = (£o(0)) 2 (£ vg ,o) 2 

:evg,o) 4 7 (fi Vii { m s , 

1 + v /l + 4(e vg> o) 4 7 ( ' R) ki| 1 + \f^ r ^v S ,o) 4 J {R) \"2\ 

— 1/2 

where £ V g,o = t vg o is the mean field VG correlation length in unit of re, and eq.(2.7) was used 
To go beyond the mean field analysis for studying G vg and enter the quantum VG critical regime 
present in T < T cr , we invoke the ordinary scaling hypothesisi'Q'S for the VG correlation function 
on the basis of the above mean field expression as follows: 

Xvg(k;^i,u, 2 ) = C9 (evgao(0)) 2 5 x (ke v g;7 {fi) (evg) z |^i|,7 (fl) (evg) z k 2 |), (2.25) 

where the correlation length £ vg is given by eq.(2.16), and the scaling function S x and a positive 
coefficient c g , as well as in eq.(2.16), may depend on Xp according to the action (2.13). 

Now, the terms corresponding to G vg in r.h.s. of eq.(2.21) will be examined. First, closely 
following the analysis§B0' used for the thermal glass transitions, let us consider the contributions 
G4 g a ^ of the diagram Fig.l (a) and of a sum of the family of Fig.l (b) to eq.(2.21) with the 

coefficient (2.22), which are expressed as 

RqG™ = ^t b I k'e~^f u (k 2 ) (-^j r 1 E Xv g (k, U , U + «) 

= ^ p c 5 (xf))/ 11 (0)Cg 2 (7 (fi) ^ g r 1 ) / k 2 Y,S' x (*; 2vr|n|( 7 ^)^ g r 1 )), (2.26) 
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= 4VTCQ1 c s (xf )) 4 (7^/T 1 ) E / <%(k; avrln^^)^/?- 1 ), (2.27) 



/k 

',2 , 



where coi = J k k 2 (f 00 (k 2 )f u (k 2 ) - (/ 01 (fc 2 )) 2 ) e" fc2 / 2 , and E^^(k; |w|) = -(d/d\Sl\) 
J2uj "Sx(k; M> l w + ^|)|n=o- Although we have assumed in eqs.(2.26) and (2.27) the "pinning lines" 
unrelated to the VG susceptibility to carry not but U p /2irr 2 3 , this simplification does not affect 
the conclusion given below, that RqG^ in the dirty limit takes a universal value independent of 
T at B vg , because of the relations (2.11) and (2.34) (see below). 

So far, no microscopic (electronic) model leading to the dissipative GL action has been specified, 
and a T = limit of the resulting GL action has been simply assumed. Although we will consider 
below a simplified BCS hamiltonian with a short-ranged repulsion and a nonmagnetic potential 
disorder in order to see how a T-insensitive GL action can be realized, it is instructive to first start 
from the dirty limit with no electron-repulsion. In this case, the expressions of coefficients U4, 7, 
/x(0), /ii(0), and C(0,0) are available in the literatureslllll) and given by 

U 4 = 87rr 3 (/3AT(0))- 1 ]T(r(2e; B)f, (2.28) 

e>0 

7 = 4vrr 2 /?- 1 E( r ( 2e ^)) 2 ' ( 2 - 29 ) 



e>0 



T 



/,(0) = In [ — ] + 47TT/T 1 ]T( r ( 2e ; °) " r ( 2e ; B ))> ( 2 - 30 ) 



e>0 



M i(o) = mo) + 47TT/T 1 E( r ( 2e ; B ) - r ( 2e ; ( 2 - 31 ) 



e>0 



and 



C(0, 0) = 2ttt(3~ 1 J2 r ( 2e ; - B ) r (2e; 35), (2.32) 

£>0 

where r(2e;i?) = (2|e|r + irl 2 B / (fto)^ 1 , and e denotes a Matsubara frequency for fermions. Note 
that the relation (2.22) is satisfied just on the H C 2{T) line where /i(0) = 0, implying that, in 
the dirty limit, the renormalized //i(0) should approach 47tt/3~ 1 X^ >0 (r(2e; B) — T(2e;3i?)) with 
decreasing B. In T < T™ f , the T-dependence of T(2e; B) is cut off by the 5-dependence, and the 
above coefficients become insensitive to T in the manner U4 — > 4r 2 r^/(iV(0) Z 4 ), 7 — > 2r(rs//) 2 = 
7 (°)(T = 0), /x(0) -> m(£/iT c2 (0)), and W (0) -> ln(3S/iJ c2 (0)). On the other hand, the pinning 
strength U p is known only in T < T™ f and given by§> 



4 



(2.33) 
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In terms of these T-independent GL coefficients, a "T = 0" critical field B vg can become well- 
defined within the dirty limit. Further, at lower temperatures than T cr which is estimated in dirty 
limit as ~ AttT^/(E f t) (oc B), we have the relation 

x p = yx 2 . (2-34) 

As a result of eq.(2.34), Xp depends only on x p . Just at B vg , the combination 7^(£ V g) z /A a s 

( R) 

well as x p and x p , becomes a constant independent of T and of material parameters according 
to eqs. (2.11) and (2.14). Therefore, G^^B = B vg ) is a universal constant divided by Rq in T < 
Min(T cr , T™ f )- Similarly, RqG^ (B = B vg ) becomes (7T/V A^)) 2 / 2 multiplied by a universal 
positive constant. 

It should be mentioned that we cannot verify directly whether the (universal) constant coi is 
positive or not, because the functional forms of /oo> fu, arid /01 are not completely known, although, 
by definition, the VG contribution G vg to the conductance must be positive. On the other hand, we 
have a unlimitedly large number of diagrams contributing to RqG vg (B = B vg ) in the same way as 
Fig. 1(b). All of them can be seen as such diagrams© that, according to the already-mentioned LLL 
diagrammatic rule, the vertex correction unrelated to the VG susceptibility is of higher order in x p 
and X4 (in T < T cr ) compared with those in Fig.l. In this sense, Fig. 1(a) is the lowest order term, 
and the diagram in Fig. 1(b) is the next lowest order term in x p . However, just at B = B vg where x p 
and X4 take constant values, the diagram of Fig. 1(a) becomes of the same order as that of Fig. 1(b) 
except the extra power in ^~ g . The same thing holds in the (formally) higher order diagrams in 
x p and »4, and hence, they also contribute to a universal value of RqG vg (B = B vg ) together with 
a sum of the family of Fig. 1(b). Since, unfortunately, we have no resummation scheme, useful at 
the critical point, for judging which of those diagrams should be adopted or may be neglected, it is 
difficult, as in the argument!! © on a T = critical conductance RqG s (T = 0,u — > 0), to estimate 
here a concrete value of RqG vg (B = B vg ). We can just conclude that, in the present dirty limit, 
i?qG vg (T — ► 0, B = -B V g) is a universal positive number. 

Of course, a universal critical G vg obtained above is not a consequence of the T = scaling 
argument ■HE!' Actually, the w-summation in eq.(2.27) was not changed above into a frequency- 
integral because 7^£ V g/3 _1 is finite in T — ► limit and at B vg , implying that G^ g \B = B vg ) is 
not a dc limit of a T = conductance!'© but a dc conductance in the quantum regime < T < T cr . 

However, it will be difficult to explain available resistivity data in terms only of the results in the 
dirty limit with no electron-repulsion. As examined in ref.19, the mean field value -B V g,o calculated 
in the dirty limit seems to increase with increasing (Eft) -1 due to the corresponding increase 
of H^ 2 {0), while the data in ref.4 have shown a trend opposite to this. Although the resulting 
fluctuation-corrected field B vg ~ £? vg) o(l — 2cr/(itEft) ) (see the sentence below eq.(2.14)) may 
decrease with increasing (Epr) -1 depending on E-pT- values, it will be difficult to understand a 
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dependence© of B vg on the film thickness d ~ kp 1 Rq / \R r E-pr) in the dirty limit with no electron- 
repulsion, where R r is the high temperature sheet resistance. Once the interplay between the 
electron-repulsion and disorder is taken into account, however, H C 2(0) and hence, B vg decrease© 
with increasing the electron-repulsion strength© Ai ~ Rt/^ttRq). Further, a nonzero Ai results 
in a failure of the equality (2.34), and thus, the G vg (T — * 0) value at a critical field is not universal 
any longer but will depend on Ai. 

We argue that the properties in the close vicinity of T = (corresponding to the region below T rep ) 
of disordered thin superconducting films have not been examined so far experimentally. According 
to the previous works, the interplay between an electron-repulsion and disorder appears in 
the GL action in two different ways: The GL coefficients C/4, n(0), fJ>i(0), C(0, 0), and U p are 
convergent in low T limit at each order of the Ai-perturbation series, and their T-dependences 
are controlled, through the denominator of Cooperons, by the |e| (Matsubara frequency) value of 
the order Z 2 /(87rrr^). Namely, their T-dependences are lost, as well as those in the dirty limit, 
below T™ f independent of Ai. In contrast, 7 at low enough T is expande in powers of 

Ailn(T/T™ f ) because this quantity is dominated by the lowest |e| values, and consequently, a In- 
dependence of 7 induced by the electron-repulsion will become remarkable rather near T = 
below© T rep ~ T™ f exp(— A^ 1 ). Since, as far as we know, R r < Rq, or equivalently 87rAi < 1, 
is satisfied in real thin films with an FSI behavior, T rep will lie much below T™ and seems to 
be inaccessibly low in real systems. Then, it is reasonable to assume the FSI behavior seen in 
real experiments to be a phenomenon in the intermediate region T rep < T < T™ f , where all GL 
coefficients are insensitive to T so that an apparent VG critical field© B* g and a nonuniversal 
constant RqG vs (B = B* g ) are well-defined. 

In §4 of ref.19, a computation result on 7 was given supporting the argument on the presence of 
the intermediate temperature region, and it was suggested that the temperature scale corresponding 
to T rep will be below 0.1 T™ f in the cases with realistic i? r -values. Unfortunately, it was difficult to 
judge whether 7 remains positive or vanishes in low T limit, i.e., in T < T rep . We simply expect here 
that 7 will significantly decrease on cooling below T rep and hence that, within the model of purely 
dissipative dynamics, the true B vg will lie at a much lower field than B* g (see Fig. 3 below).© What 
we wish to emphasize is that, as far as 87rAi < 1 is satisfied in real systems, the T-dependence of a 
microscopic origin is remarkable rather at extremely low temperatures below T rep , and hence that 
any attempt, such as ref.6, to explain the FSI behaviors by assuming a T = bosonic model is not 
justified. 

§3. Description of 2D Resistive Behavior near B* g 

In this section, we discuss the resistivity curves around B* g and in T < T™ on the basis of the 
results in §2 and give some results relevant to comparing with experimental data. Again, the GL 
coefficients will be assumed for a moment to be insensitive to T so that a critical field B vg may be 
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well-defined. 

First, let us start with the resistive behaviors far above B vg [B > B vg ) where the Gaussian 
approximation, illustrated as eq.(2.24), for the VG fluctuation may be used. Assuming the low 
frequency behavior to be essential even in the Gaussian region, we will keep, for simplicity, only 
0(|u;j|) (j = 1, 2) terms in eq.(2.24), and this Xvg,o will be substituted into the first line of eq.(2.27). 
By arranging the w-summation and performing the |f2|-derivative, we obtain 



xln| 1 + L 7 ( R )g 4 ' +ln ( 1 + c c £v g ,o)^— (6o(w) + 0o(O)) 



-(SoH) 5 



(3.1) 



which obviously vanishes in T — > (/3 -1 7^£vg,o — * *-*) limit. Similarly, one can verify that G^gQ 
also vanishes at T = if taking account of the frequency dependence of Q\. The additional 
ln^gQ dependences arise from an upper cutoff {c c rs)~ l (with a constant c c of order unity) of 
the |k|-integral. Judging from a similar situation one encounters in deriving 2D AL fluctuation 
conductance in T — > limit and at B = 0, we believe that this divergence is specific to the present 
direct frequency-summation and may be avoided by the standard analytic continuation which we 
have not tried. However, this technical issue does not affect our conclusion that, as well as Gfi, 
G vg (B > B vg ) vanishes in T — > limit, because, as shown in ref.14, each term of perturbation series 
of fluctuation conductivity examined within a quantum GL action vanishes in T — * irrespective of 
the presence or absence of a pinning-disorder term in the action. In B ^> B vg , such a perturbation 
series should become more convergent with approaching T = 0, and hence, we can conclude that, 
as well as each term of the perturbation series, the resummation results of the perturbation series, 
i.e., Gfi and G vg themselves in B > B vg , also vanish in T — > limit. 

On the other hand, in B < B vg and out of the quantum critical regime defined by eq.(2.15), 
£ V g grows with decreasing T or B and hence, the VG fluctuation becomes "classical" even at 
low temperatures below T cr in which the pair-field fluctuation, with higher energy than the VG 
fluctuation, is of a quantum character. In fact, as eq.(3.1) suggests, one may expect the 2D classical 
behavior!) G vg ~ /3 -1 7^£vg within the present analysis. However, the classical (i.e., thermal) 
2D VG transition is washed out, e.g., by a nonperturbative effect such as the free vacancies or 
interstitials in the vortex solid© Therefore, the region in which the classical scaling behavior 
G v g ~ is visible may be quite narrow. 

According to eq.(2.27), the scaling behavior of RqG vg 

R Q G vg =u(c u ^-^(^\ ) (3.2) 
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is expected in the quantum VG critical regime denned by the relation (2.15). Here, U(x = 0) is a 
positive and nonuniversal constant according to the result in §2, and 



Further, if describing the "Gaussian" region, discussed above, in terms of eq.(3.2), the limiting 
behaviors lA{x —* +oo) — > and U(x — > —1) — > (x + l)~ zv have to be satisfied. Therefore, the 
resulting {RqG vg )~ l v.s. T curves around B vg in B-T plane are, just like the original speculation© 
based on some B = results, similar to that of a " renormalization-group flow" near a fixed point. 
Further, we note that, if zv > 1, the 7-dependence of c u will be opposite to a naive expectation 
on a strength of quantum critical fluctuation. In the case with a quantum normal-Meissner (i.e., 
B = 0) transition, for instance, one would expect the width of the quantum critical regime at a 
fixed T to scale like "f l l zv and hence to become narrower for a stronger quantum fluctuation. In 
contrast, in the present B > problem, the quantum VG critical region becomes wider as the 
pair-field fluctuation is enhanced. Details of the low T regions in B-T phase diagram of thin films 
are sketched in Fig. 2. 

As shown previously]^ RqGji calculated in the pinning- free (U p = 0) case shows, around a field 
above B m (0), "fan-shaped" resistivity curves similar to but more moderate than that of G vg near 
B vg . This follows from the facts that, when 7 and hence the melting line B m (T) in the quantum 
regime are insensitive to T, a "flat" Gj t curve of the order of Rq is realized at a field above B m and 
that, in higher fields, Gj^ shows an insulating behavior reflecting Gfi(T = 0) = 0, while it yields, 
in lower fields closer to B m , the classical vortex flow behavior which is insensitive to T at such low 
temperatures (see, for instance, Fig. 4 of ref.14). Although an inclusion of the vortex pinning effect 
will slightly change this behavior of Gfi particularly close to B m , it is generally questionable to 
interpret real resistance curves by neglecting the presence of G fi and identifying the resistance only 
with G~ g because, as mentioned above, Gfi can have a magnitude of the same order as G vg near 
B m . The fluctuation corrections to G n , 5G S , will be taken into account later, and for a moment we 
identify the total conductance with the bosonic contribution Gfi + G vg plus the normal contribution 
G n . 

Now, we discuss typical resistance data under the assumptions that T rep is inaccessibly low (see 
the end part of §2) and that an apparent critical field B* g , well-defined in T > T rep , will lie above 
B m . Due to the former assumption, all the GL coefficients and hence, B* g are insensitive to T 
in T rep < T < T™ { . One will see soon that the latter assumption has already been verifiedBB 1 in 
some data. First, let us start with the cases with weak disorder (small R r ), in which the quantum 
fluctuation is also weak simultaneously and B* g will lie near but above B m (see below). In this case, 
the contributions of G n and/or Gfi occupy a large weight of the total conductance near B* g , and 
there may be no clear indication of quantum correction!!) (~ — Ailn(l/Tr)) to the familiar residual 




(3.3) 
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behavior (insensitive to T) of G n at accessible temperatures above T rep . Since T™ f r < Bl 2 /cj) < 1 
in fields of our interest, the negligible quantum correction to G n implies that the neglect of the 
region below T rep is justified. Further, since Gfi just above B m in the present case is, as mentioned 
in the last paragraph, insensitive to T, the background contribution G n + Gfi is expected to be 
independent of T and to weakly depend on B. On the other hand, according to eq.(3.3), the 
quantum VG critical region at a fixed T seems to be narrower at weaker disorder: As a rough 
estimation in this weak disorder case, if any Ai-dependence is neglected in eq.(3.3), c u certainly 
decreases with increasing disorder as far as zv > 1, as seen in most of data. By taking account 
of these contributions to the conductance altogether, the FSI behavior (3.2) is expected to be 
visible in a narrow VG critical field range around B* g . We believe that this will be an appropriate 
explanation, for instance, to the data in MoSi filmsi) where the total resistance value near B* g was 
remarkably suppressed due to a large weight of Gfi + G n insensitive to T. In fact, in ref.4, the field 
Bq below which the thermal activation barrier!!) ~ ml?" 1 , arising from the motions of vacancies 
and interstitials in a (short-range ordered) vortex solid, remains nonvanishing lies just below the 
critical field corresponding to B* g . Since Bq should be essentially the same as B m , it justifies the 
above assumption that B* g lies just above B m . 

It is possible that, in systems with still weaker disorder, even the quantum fluctuation behavior 
of Gfi is not seen. In such a case with inaccessibly low T cr , the flat behavior of G vg should not be 
seen consistently. This situation corresponds to the data in ref.33 where the only quantum behavior 
was seen in a vortex flow behavior, which is itself the classical behavior of Gn, but suggestive of 
a quantum tunnelling effect. Clearly, the metallic behavior below H C 2 there!!* is not a reflection 
of a T = phase diagram: When the sample disorder is weaker, low enough, one needs 

to enter a lower temperature region to find quantum VG behaviors. Theoretically, an intermediate 
vortex liquid phase with finite resistance much smaller than G~ l is absent^ in 2D homogeneously 
disordered films at T = (see also §4). 

Returning again to systems with visible quantum critical behaviors of G s , in turn, let us consider 
situations with stronger disorder (larger R r ). In this case, a decrease of Gfi and G n on cooling 
will affect the total conductance: The total conductance at B* g will decrease on cooling until a 
constant limiting value G ~ G vg (B* g ) is reached. Such a behavior that the resistance just at 
B* g is not flat but increases on cooling until the lowest (accessible) temperature is reached has 
been observed in various materials.!'© If such data are explained according to this scenario, the 
vanishing contribution G — G vg in T — > will not be accompanied by the scaling behavior (2.15). 
However, an insulating behavior at B* g may appear even if G ~ G vg : As illustrated by eq.(2.26), a 
sub-leading term of G vg {B; g ) will behave like ~ (f vg (B = B* g ))- 2 oc T 2 I Z . If the sum of such terms 
is positive just like eq.(2.26) itself, this also becomes an origin of an insulating behavior at B* g . 
We note that, in this case, the sub-leading terms in the quantum VG critical region are multiplied 
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by a critical scaling function like U{x) in eq.(3.2) and hence that they can be discriminated from 
an insulating behavior of the noncritical term G — G vg mentioned earlier. 

In ref.7, nonmonotonic resistance curves (insulating at intermediate temperatures but supercon- 
ducting on further cooling) were found even just below B* g in Bi/Sb films. According to our 
theory, this nonmonotonic behavior below B* g is also understood as arising from the sum of the 
above-mentioned insulating Gj^ above B m (0) and the superconducting G~ g below B* g . This non- 
monotonic behavior is visible because the window in B m (0) < B < B* g is moderately wide in 
contrast to that in MoSi caseS Actually the presence of a wide region B$ < B < B* g , in which 
the resistance does not yield an activated behavior indicative of a vacancy- or interstitial-creep in a 
vortex solid, was pointed out in ref.7. As also mentioned in §3 of ref.19, this window should broaden 
with increasing disorder. However, our theoretical result shows that it is invalid to, based only on 
the nonactivated resistive behavior, identify^ the window in B m (0) (or Bq) < B < B* g with a pu- 
tative quantum liquid phase with low but nonzero resistance at T = 0. Of course, in a system with 
still stronger disorder or at low enough T (but above T rep ), the contribution G n + Gfi is already neg- 
ligible, and the leading term of G vg (B = B* ), i.e., a (nonuniversal) constant G ~ G vg (B = B* ), 
should be observed as the total conductance, as actually seen in the state 1 of ref.3 where a critical 
resistance value is much larger than R r . On the other hand, the apparent critical resistance values 
in ref.7 seem to correlate to the corresponding i? r -values, implying that the samples in ref.7 have 
intermediate strengths of disorder. 

We emphasize again that, in the present theory, the intervening metallic behavior at a field 
is regarded not as a reflection of the true T = phase diagram but as a phenomenon in the 
intermediate (accessible) temperature range T TC p <C T <C . This interpretation is never artifitial. 
In fact, a 2D FSI behavior was also observed in 2D-like but bulk (underdoped) YBCO,& in which 
the true T = critical behavior must be of a 3D VG type (see §4), and the system at nonzero 
temperatures behaves as if it have a 2D VG transition at T = (This situation is comparable 
with a dimensional crossover just above a thermal transition. Note that, in the present case, 
Cvg(B = B* g ) diverges in T —* limit). Similarly to this, in the present case the temperature 
variation of resistance curves in B > B vg (see Fig. 2) should become insulating at inaccessibly low 
temperatures below T rep , although the resistance curves in B vg < B <C B* g may decrease on cooling 
in T > T re p, more or less, as a result of classical (thermal) VG fluctuation. 

Now, the fermionic fluctuation term 5G S of the conductance will be considered. Since discussing 
the normal part G n in details is beyond the scope of this paper, let us assume here G n {T > 0) to 
show an almost metallic behavior. Then, the dynamics of the pair-field should be dominated by the 
dissipative term, and the fermionic fluctuation conductance 5G S , consisting of the Maki-Thompson 
terms and DOS terms, remains nonzero in low T limit and was previously regarded as a correction 
to G n without examining in details (see §7 in ref.14). According to a recent systematic study in 
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ref.25, SG S at T = is negative and given, on the Gaussian level, by 

which is valid above the mean field -ff C 2(0) (/^(O) ~ — 1 + B/H^ 2 {0) > 0). Note that, in the dirty- 
limit, the r.h.s. of eq.(3.4) is independent of material parameters. If the renormalization of LLL 
fluctuation is performed in terms of eq.(2.5) in deriving the corresponding one to eq.(3.4) applicable 
in B < H C 2(0), one finds eq.(3.4) is replaced by 

^ g .H> (T = ),- ^""(;;°^ K0)|, (3.5) 

where effects of pinning disorder on the fluctuation renormalization were neglected for simplicity. 
By using eqs.(2.28) and (2.29), the r.h.s. of eq.(3.5) becomes of the order of — Epr\[J,(Q)\. Further, 
the increase of \SG S \ on cooling© remains valid after the fluctuation renormalization is performed. 
Since |^(0)| < 1 within the GL theory, while RqG u ~ EpT in 2D and with no quantum correction, 
the fermionic conductance G n + 5G S can significantly reduce with decreasing B below H C 2(0) 
due to the |/x(0) [-dependence in eq.(3.5). This seems to explain the negative magnetoresistance 
(MR) in MoSi thin films, which was not visible in 3D-like films and in highly disordered films 
nonsuperconducting even at B = 0.0* Further, a more remarkable negative MR had been also 
found in InO thin filmsSS showing the FSI behavior. In ref.34, the origin of this negative MR 
had been ascribed to a localization of bosons (i.e., pairs) which is also of a superconducting origin 
but, in contrast to our idea, does not seem to be supported through a microscopic calculation. We 
also note here that the data in ref.34 have suggested a metallic resistance in much higher fields 
than B* g and even at low enough T. In our notation, this corresponds to a nonvanishing G n and 
is consistent with our assumption that, in contrast to ref.6, the FSI behavior should be explained 
based on a dissipative dynamics of the pair-field. 

Our scenario on the FSI behaviors is conventional and precludes a possibility of an intermediate 
metallic vortex phase at T = 0. An argument favoring such an intermediate phase is based on 
the data suggestive of a quantum tunneling behavior in a lnG -1 v.s. 1/T plot!) and also on 
a computation of fluctuation coductivity applicable only to B = and based on a neglect of 
dissipative dynamics© As verified in ref.4, however, such a metallic behavior over a broad field 
range was also seen in the case irrelevant to the vortex states, i.e., the case in a field parallel to the 
surface of thin film samples and with a current parallel to the field. This finding suggests that the 
the quantum tunneling behavior!) may not be due to an intrinsic origin. Actually, an intermediate 
metallic behavior tends to be seen in rather weakly disordered films with R r < 1 kfi. As mentioned 
earlier, one needs to enter lower temperatures in a weaker disorder case in order to search a true low 
T behavior. To the best of our knowledge, there is no clear evidence of an intermediate metallic 
behavior in a stronger disorder case. Actually, an intervening metallic phase was argued based 
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only on models for the granular case in zero field which are inapplicable to the nonzero field case 
of homogeneous materials (see §1). Further, the argument® favoring G(B = B c ) ~ G n certainly 
contradicts the -R r -dependence of critical resistance values measured in ref.3. 

§4. Discussions and Extension 

In §3, we have discussed various resistive data suggestive of the FSI transition at T = by taking 
account of three different superconducting terms Gfi, G vg , and 5G S of conductance. Those data 
have been explained by noticing that the experimentally accessible temperatures in real thin films 
with R r < Rq will lie not in the vicinity of T = but within an intermediate temperature range 
defined by us. The resistance data in various samples are apparently incompatible with one another 
and are not explained comprehensively within a scenario taking account onlyi) of G vg without a 
detailed calculation. As in the thermal casef© a couple of contributions to conductance are 
necessary in order to explain avalable data in a unified manner. Although we have explained, by 
focusing primarily on systems with relatively weak disorder, why the apparent critical resistance 
value tends to increase with increasing R r , it has not been clarified whether the G vg (B = £>* g )-value 
itself decreases or not with increasing R r . Further, attention was not paid much to the values of 
exponents z and v. Although they are usually assumed to take universal values, this assumption is 
not necessarily valid for a quantum transition in random systems. These issues are left for future 
studies. 

The extension of the present theory to 3D case will be explained here since, at a glance, recent 
dat;S in thick (3D-like) MoSi films would seem to contradict the present theory. Similarly to the 
2D case, the conductance G± for a current perpendicular to B can be seen as consisting of the four 
terms; 

G± = G nj ± + Gji t ± + 5G S) ± + G vg)J _, (4.1) 

where the notation follows that in 2D case. Although the conductance defined in the applied field 
direction is not discussed here, our final conclusion remains valid for the parallel conductance. First, 
let us list the behaviors of each component in eq.(4.1) near T = 0. In 3D, the normal part G n ± is 
safely assumed to be metallic, and the fluctuation correction 5G S ,± is nondivergent in contrast to 
2D case and merely a small correction, insensitive to B and T, to G n , i.e., 5G S ,± ~ G ny ±_ 0(1/ Eft). 
According to ref.14, the AL fluctuation term Gfi ± (except G vgj j_), as in 2D case, approaches zero in 
T — > limit above any transition field (see below). Further, it is easily found as a trivial extension 
of eq.(2.27) that the VG contribution G vg) ± has the form of a scaling function multiplied by ^g 1 , 
and thus, 

G vg (B vg ) ~ C/(r, B vg ) ~ T 1 /*, (4.2) 

which implies G vgt ±(T — ► 0, -B vg ) = 0. Thus, above any transition field and at T = 0, the total 
conductance is essentially equivalent to G n . How about the B-D phase diagram (Fig.3) at T = 0, 
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where D measures the strength of pinning disorder ? In this case with dissipative dynamics, 
the quantum 3D GL model is equivalent to a classical 5D one, and the dimensionality of LLL 
fluctuation is three. Thus, the ordinary critical line H* 2 (T = 0), signalling the onset of the ordinary 
long-ranged phase coherence, lies in nonzero fields in this 3D case, while the corresponding field 
in 2D is zero. Just as in the B = thermal transition in a bulk superconductor, the ordinary 
critical point H* 2 (0) (the left vertical line in Fig. 3) is lowered from the mean field H c2 (0) (the 
right vertical line in Fig. 3) due to the 3D quantum fluctuation. When the pinning disorder is 
absent (i.e., D = 0), as mentioned previously in §5 of ref.14, a first order vortex solidification 
transition at T = and in 3D should occur above H* 2 (T = 0). Its position, the end point of 
5 m (0)-line, was indicated as an open circle in Fig. 3. On the other hand, as found in ref.37, the 
ordinary superconducting transition in classical 5D case with nonzero D is expected to be of second 
order. This fact may be useful in understanding the present situation due to the above-mentioned 
correspondence between the classical 5D and the quantum-dissipated 3D cases. Since the first 
order solidification at the open circle is also accompanied by the ordinary superconducting ordering, 
however, both transitions must connect with each other in B-D diagram, as described by the chain 
line. Then, it is straightforwardly concluded that the resistance for any current direction must 
be zero at least at and below Max(f? m (0), H* 2 (0)): According to the definition of VG ordering 
(see eq.(2.10)), the presence of the ordinary phase coherence inevitably implies the presence of 
VG ordering, although a VG ordering can generally occur© with no ordinary phase coherence. 
Namely, if an estimated -B vg (D)-line lies below H* 2 (0) for some .D-values, then the resistance for 
any current direction vanishes not at B vg but already at H* 2 (0). Thus, since there is no region 
with finite resistance below H* 2 (0), no metallic intermediate phase with much smaller resistance is 
possible in 3D and at T = 0. Although, at large D, B vg may lie above H* 2 (0), it is obvious that 
this conclusion is not affected. The resistance curves apparently residual on cooling, observed in 
ref.36, are the reflection not of a true T = metallic vortex phase but, just like the data in ref.33, 
merely of a pinning-induced enhancement of the vortex flow conductance in the vortex liquid region 
becoming narrower in T —* 0. 
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(a) (b) 



Figure Caption 

Fig.l 

Examples of diagrams contributing to G vg . The solid curves denote the LLL fluctuation propaga- 
tors, the chain curves are NLL propagators, the dashed lines are the pinning lines, and the hatched 
rectangle denotes Xvg- 

Fig.2 

Details of 2D phase diagram near zero temperature. Here H C 2(0) is the mean field upper critical 
field, B m (0) the T = vortex-solidification field estimated by neglecting the presence of the critical 
field B vg , the chain line denotes T rep , and the lower (upper) dashed straight line denotes T™ f (Z?) 
(Tcr(B)). The quantum VG (QVG) critical regime is the T-dependent field range defined by 
eq.(2.15) around the apparent VG critical field B* g . Details of the true zero temperature limit 
below T rep are not drawn. 

Fig.3 

Conjectured 3D B v.s. D (pinning disorder) phase diagram at T = 0. The chain curve B m (0) 
and the open circle imply the first order transition accompanied by a melting of a vortex solid, 
the left vertical line is the ordinary superconducting transition line H* 2 (0), the right vertical line 
the mean field H C 2 line, and the curved portion in large D of the solid line indicates B vg occurring 
above H* 2 (0). The solid lines imply second order superconducting transition. As explained in the 
text, the VG susceptibility is divergent everywhere below Bm = Max(ff* 2 (0), B vg , B m (0)) and 
hence, the resistance R± is zero in B < Bm, while R± ~ G^ 1 in B > Bm- 
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